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Abstract 

A one-dimensional equation describing the distribution of the effective vector potential A{y) across a film 
width, which holds for thin (d < A) and thick (d > A) films alike, is derived based on the analysis of a 2D 
Maxwell-Londons equation for superconducting films in a perpendicular magnetic field. The validity of this 
equation for a finite-thickness film is verified by a numerical analysis. An approximation dependence A{y), 
finite (with all of its derivatives) across the entire film width, is found for films, being in the Meissner state. 
The flux-entry field is evaluated for a film of arbitrary thickness. An approximation expression is obtained 
for the distribution of the sheet current density in the mixed state of a pin-free superconducting film with 
an edge barrier. The latter approximation allows to estimate magnetic field concentration factor at the film 
edge as a function of external magnetic field and geometrical parameters of the sample. 
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1 Introduction 

Of active interest recently is theoretical ||, ||, ^, |6[ 0, ||, ^ and experimental 10 investigation of 
mixed static and dynamic states in superconducting films in a perpendicular magnetic field. Theoretical cal- 
culations of various magnetic characteristics of films in such a geometry, using the microscopic theory or the 
Ginzburg-Landau equations does not seem possible because of the mathematical intricacies involved in their 
solution. Therefore, in practical calculations the following two approaches are mostly employed. In the first 
method the Maxwell equation is analysed jointly with the London equation, which yields an integral-differential 
equation [0, ^, 0, § for distribution of the current density integrated over the film thickness d. This 
equation was derived on the assumption that the film is thin (d <C A, A is the London penetration depth), which 
naturally limits its applicability range. 

The other approach is based on the theory of complex variable functions (TCVF) |j, used in transfor- 
mation of integral equations In this case the phenomenological dependences B(H) and J(E) are employed 
as an additional condition instead of the London equation. It should be noted that neither of these methods 
has ever focussed on investigating the effects related to finiteness of the London penetration depth A. Besides, 
the equality of the results obtained through solving of the integral equation for thin films and by the use of the 
TCVF methods for thick films leads us to believe that for films of arbitrary thickness there should exist one 
equation describing distribution of the current density and the magnetic field across the film width. 

The present paper deals with a study on the distribution of the magnetic field and current density in the 
Meissner and mixed states for films placed in a perpendicular magnetic field. It is shown that for the finite- 
thickness films the Maxwell-London equation ^ describes distribution of the vector potential (current density) 
averaged over the film thickness d, provided the latter, is much smaller than the film width W: d/W <C 1. For 
thin films, d ^ A, this equation is valid practically over the entire width of the film, while in the case c? A it 
holds everywhere in the film, except for the areas near the edges, which measure as W/2 < \Y\ < W/2 — d/2. 
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The paper is organized as follows. Section 2 describes derivation of an one-dimensional equaton for dis- 
tribution of the film-thickness-averaged vector potential across a sample width, based on the analysis of 2D 
Maxwell-London equation. In subsections 2.1 and 2.2 the 2D and ID equations arc numerically studied and 
compared for thin and thick films respectively. An approximation dependence for A{y) is obtained through 
numerical solution of these equations. The distribution of the vector potential (or the local current density) and 
of the local magncti field over a superconductor cross-section is described. Section 3 deals with estimation of 
the field for the first vortex entry into thin and thick superconducting films (subsection 3.1 and 3.2 respectivly) . 
In subsection 3.3 we discuss the influence of surface defects and a layered structure of superconductors on a 
barrier suppression field valu{\ Sec;tion 4 is concerned with the structure of a mixed state in thin and thick 
superconducting films of an arbitrary width in the absence of bulk pinning. An approximation formula is found 
for the integral (over thickness) current density, which is then used as a basis for constructing the magnetization 
curves for these superconductors. Section 5 sums up the main results obtained in this work. 



2 The structure of the Meissner state 

Assume an infinite (in the X direction) superconducting film of width W and thickness d in a perpendicular 
magnetic field; the geometry is shown in Fig.l. Let us first consider the Meissner state. The Maxwell equation 
has the form (in a gauge V • A = 0) 

47r 

AA = -^j, (1) 

where A is the 2D Laplacian operator. It also follows from the symmetry of the problem that only the 
X components of the vector potential A = {Ax,0,0) and of the current density j = {jx,0,0) are not zero. 
The boundary conditions are ^^\y^±oc = — -ffoo, ^|^|z^±oo = 0, where is the field far from film 
H = (0, 0,i?oo)- It should be noted that by the magnetic field H in a film here we imply a microscopic field 
averaged over scales much larger than the atomic one but much smaller than A. 

In this case it is convenient to change over from the differential equation (1) to its integral analog. Using 
the Green function of the 2D Laplacian operator, we rewrite Eq.(l) as 

A,{Y,Z)=Al{Y)-- / {\n\K--R'\+C)j,{Y',Z')dY'dZ', (2) 

C J -w/2 J-d/2 

where C is the constant generic for the 2D Green function, A'^ (Y) is the vector potential of an external field: 
A^{Y) = —HocY and the origin of coordinates is chosen in the film centre. 

Employing London equation j = — A/47rA^ and introducing dimensionless coordinates y = 2Y/W, z = 2Z/d, 
Eq.(2) reads 



A4y, z) = Al{y) + ^ In (^y - y'f + ' (z - z'f^ A.{y', z')dy'dz' 

Wd ~ 

+ Y^C' Ax{y', z')dy'dz', (3) 

where C = C + 2\n{W/2). The latter integral in (3) is directly proportional to the total current. In a 
magnetic field (without a transport current) the total current is equal to zero due to the symmetry, so the last 
term in the right-hand side of Eq. (3) vanishes. We now average Eq.(3) over the film thickness, which yields 
the following expression (hereafter by A{y) we imply 0.5 Ax{y, z)dz) 

A{y) = -H^yW/2 + ^ £^ In \y - y'\A{y')dy' 

+ 3^^ /_i /i /i ^""^^ + ^'/^^ " y'^'^^-^y'^ z')dy'dzdz', (4) 
where the integral kernel of Eq.(3) is written in the form 
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and the designation b = {z — z')d/W <C 1 is introduced (as in this case of superconducting films d/W ^ 1). 

The function ln(l + b"^ / {y — y'Y) is not zero only in a small region around point y': \y — y'\ < \b\. In this 
case integration over y' in the second integral of Eq.(4) can be done only over this small region. For the same 
reason we can expand the function Ax{y' , z') into the Taylor series in terms of y' near point y: 

A.{y',z') = A,{y,z') + ^^^^^{y' - y) + ■ (5) 

Note that expansion (5) (in the above specified limit) is valid for thin [d < A) films over the entire sample 
width. In thick {d > A) films the validity of this expansion is violated in the near-edge regions with dimensions 
of order d/W. More details on the applicability of Eq.(5) are provided in the end of this Section. 

After the series-expansion of function Ax{y', z') (5) it is now possible to calculate the last term of Eq.(4) 



iT l\ ^ (^) ' " "'^''^^ " ^'^') (^^^^' "'^ ^ ^^fy"'\ y' -y) + ...) dy'dzdz', (6) 

(note, that the upper (lower) limit of the integration over y' will change by 1(— 1), when point y becomes close 
to the film edge, i.e., when 1 — |6| < || < 1; this takes place because the integration in (4) is carried out only 
over a sample volume) and to show that (6) is equal to zero in a wide parameters range. Indeed, integration 
of Eq.(6), first in terms of y' and then in z and z', provides a direct evidence (bearing in mind that function 
Ax{y', z') is even in z') that integral (6) is zero for all values of y satisfying the inequality |y| < 1 — \b\. In the 
region 1 < |y| < 1 — \b\ integral (6) leads to appearance of nonzero terms that for thin films arc small due to the 
presence of the corrections of (d/A)" (n > 1) type. They have to be taken into account when we are interested 
in, for example, the distribution of the derivative dA/dy near the film edges (since disregard for these terms will 
cause a logarithmic divergence of the first derivative). For thick films, allowance for the near-edge regions of a 
superconductor in integral (6) cannot largely affect the A{y) distribution off the film edges because of smallness 
of|6|. 

Thus, the 2D equation (3) is reduced to a ID equation for the film-thickness-averaged vector potential Aiif), 
which is valid in the region \y\ < 1 — |6| 

A{y) = -H^yW/2 + ^ £^ In \y - y'\A{y')dy' (7) 

In the following sections the results of a numerical solution of Eqs. (7) and (3) are provided for thin and 
thick films. 



2.1 Thin films {d < A) 

It turns out that in films satisfying the condition d/X< 1/4 the difference between the solutions of Eqs (7) and 
(3) (averaged over thickness) is about 1% far from film edge and less than 4% in a narrow near-edge region. 
An appreciable error in the near-edge region (which is slightly growing towards the film edge with a larger 
numerical step) depends on the presence of the small corrections in (7) that were ignored. 

Besides, the vector potential in this case is practically independent of z (but not the derivative dAx{y, z)/dz). 
Thus, at d = A the relation Axiy^l)/ A^rXy.Q) « 1.07, i.e., the difference is about 7%. With a lower d/\ ratio 
the relation Ax{y, l)/Ax{y,Q) tends to unit. 

We have derived an asymptotic expression for the vector potential distribution, satisfying Eq.(7) (and, hence, 
(3) averaged over z) with a sufficently high accuracy (not less than 3% at a film edge and far from edge region, 
and not less than 6% in the near-edge region, see Fig. 2): 

A{y) = - /.f^r,^,, . (8) 

where Xeff. = A^/d, /3 ~ 2Ae///7rM^ 4- 4(Ae///VF)^, and the dependence a{W/\eff) is shown in Fig.3 
together with the approximation (9) 
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At W < Ae/f the dependence ^(y) becomes almost linear. Formula (8) with a = 0.25, /3 = was derived 
analytically in m, by solving Eq.(7) (to be more exact, a simplified version of Eq.(7) in which the left-hand part 
is omitted, which corresponds to the condition WXefj >• 1). 

The resultant dependence A{y) allows to calculate the field concentration parameter 7 = H'^f^'^/Hoo (where 
H ^ ^ is the edge field averaged over superconductor thickness) for films of such type: 



a 



Hoo WVP V V 13 

At W Xeff Eq.(lO) transforms into 



(10) 




where the coefficient preceding \/W/Xeff is a quantity of order unit. The difference between Eq.(lO) and 
the numerically obtained expression for 7 may reach 30%: for wide films, {W ^ ^eff), Eq.(lO) yields an 
overestimated result, see insert in Fig. 4. This takes place because, unlike the A{y) function itself, the first 
derivative (8) with respect to y (magnetic field) adequately satisfies the numerical solution everywhere except 
for the narrow region near a film edge (see Fig. 4). Due to the logarithmic divergence of the magnetic field at 

a film edge, which follows from the solution of Eq.(7), the approximation expression for ^ was compared 
with the numerical solution of Eq.(3). Fig. 4 also shows the interpolation function 1 + 0.66-y/M^/Ae// for the 
numerical analysis data (see solid line in the insert). Thus, the difference between (10) and the numerical result 
in the wide film limit W/Xeff ^ 1 comes to about 17%. 



2.2 Strip of finite thickness {d ^ A) 

A comparative numerical analysis of the solutions to Eqs (3) integrated over a superconductor thickness and 
(7) was also carried out for the case A <C d ^ VK. It was found out that the solutions coincide (to the accuracy 
of about 3%) in the region \y\ < 1 — d/W and ( quite surprisingly) in points \y\ = 1. In the near-edge region 
we observed a difference in the solutions of Eqs. (3) and (7). 

An approximation equation has also been derived for the dependence A{y). It turned out to be exactly the 
same as the dependence (8) with the selection a ~ 0.25, f3 « 0.64A^/dW^. One can easily see that the obtained 
values for a and P practically coincide, to a numerical error, with those obtained for thin films at }? /Wd <C 1. 
Fig. 5 shows the dependence A{y) derived from the solution to equation (3), and also Eq.(8). It is seen from the 
latter that maximum deviation of the solution to Eq.(3) from (8) (and, hence, from the solution to (7)) occurs 
in the region y > \ — d/W, but in point y ~ I, however, both solutions coincide again. Thus, expression (8) 
provides an adequate description of the A(y) distribution (the difference from the numerical solution of Eq. (3) 
does not exceed 3%) in the region < 1 — d/W and in points \y\ — 1. This confirms the above statement that 
Eq.(7) describes well the distribution of the A{y) dependence in a film depth. Besides, Eq.(7) is apparently 
valid immediately at the edge of a superconductor as well, which, in our opinion, is quite surprising fact. 

Fig. 6 shows the distribution of the current density over a film cross-section and the distribution of the 

absolute value of local magnetic field (|_ff| = ^Jll^~\~ll^) both inside and outside a film with dimensions 

W = lOOA, d = lOA. As seen from Fig. 6a, c, the magnetic field reaches its maximum at the corners (side edges) 
of a superconductor. At the equator {y — l,z = 0) the magnetic field is less intensive, but not appreciably 
smaller than the field at the corners (for comparison, at the corners \H\ « 4.1_ffoo, while in the middle of a side 
face ~ 2.%Hoo for the given parameters of film). It is easily seen that towards the film interior (along the 
y-axis) the magnetic field is practically uniform through the entire sample thickness, except for the near-surface 
areas with the thickness of order A where magnetic lines abruptly change direction. A similar behaviour is 
demonstrated by a current density (see figs. 6 b,d). It is proved numerically that both the current density and 
the magnetic field fall off towards a sample centre by a law similar to the exponential one, not only from the 
side faces but also from the top and bottom ones. So, in thick films A ^ d <C screening currents fiow only 
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in the near-surface layers of thickness about A. On the same scale there is a decrease of a local magnetic field 
in superconducting samples of such type. 

The numerical solution of Eq.(3) also provides a possibility to determine the field at a film edge. Unfor- 
tunately, unlike with thin films, Eq.(8) differs largely from the numerical result for the near-edge region (this 
discrepancy may reach 30%). Therefore, the use of (8) in calculations of a thickness-averaged z-component of 
the magnetic field at a film edge certainly leads to a considerable error. 

In Fig. 7 the obtained numerical dependences H^f^^/Ho^ (the z-component of magnetic field, averaged over 
thickness) and Hz{l,Q) / Hoo (the z-component of magnetic field on equator) on \JW j d are presented. It is clearly 
seen that with a good accuracy the dependences are linear even for the W/d values which are close to unit. 

Note that the coefficient of proportionality between h1 /Hoc \fwjd is equal to ~ 1.03, which practically is 
the same as its estimate (~ 1) found in Generally speaking, the value of this coefficient depends on a film 

thickness (or, rather, the ratio d/A). Thus, the insert in Fig. 7 illustrates relationships H'^^^^ / Hao, Hz{l, 0)/i?oo 
for various values of film thickness and shape parameter W/d = 5. It is seen that only the quantity Hz{l,0) is 
practically independent of the ratio d/X. The strongest dependence on film thickness is exhibited by Hz{l, 1) 
(it grows with the increase of the ratio d/X). This results in a slight increase of the average field h'^^^'^ with 

a growth of film thickness (given the same W/d ratio). However, for very thick films, d ^ X^ H'f^^/Hao is 
supposed to practically cease to be dependent on d/X. Indeed, in the limit of interest the equator field is 
independent of d/X, while the corner field Hz{l, 1) increases as y^W/X (which was derived from the expression 
(11c, d) given below)). The sharpest variation of the magnetic field intensity occurs at a distance of order A 

from the top/bottom surfaces. Correspondingly, the contribution of those regions in H'^ / Hoo will be of order 
XyJW/X/d — ^JW/d^JXjd and will become negligible with the increase of the ratio d/X. 

Besides the approximation expression for A(\j), we have found numerical estimates for the vector potential 
in points (1, 1) (on a corner); (1,0) (on the equator), and also the distribution of the vector potential (current 
density) over the upper/lower surfaces: 

?W 

AM,0)-~H\l—X, (11a) 
^,(l,l)~A,(M//2,0) ( 1 + -^^ ) , (116) 




where 



A,(y, ±1) = , -^'JJ-o-^ (^^ ) 



(H-(d/A)2/27r)' " nW{l + {d/X)/ 

It is easily seen that at d ^ A the vector potential in points (±1, ±1) will largely exceed its value in points 
(±1,0). Using expression (11c) which is similar to (8) with renormalized parameters a, /3, we can find the 
points (lines) on the upper and lower surfaces, where the vector potential will coincide in absolute value with 
that on the equator. Simple calculations show that these lines are at a distance ~ d/7r from the side surfaces of 
the strip having A ^ d ^ 

These results allow to confirm the assumption (see subsection 2.1) on the possibility of expanding Ax{y' , z') 
in the limits {y — b,y + b). Indeed, in the case of thin films the vector potential (or current density in a mixed 
state; see Sec. 4) varies on scales much larger than A and, hence, than d (see (8)). For thick films, the vector 
potential (current density) far from edges is finite only in the surface layers of thickness of order A. At the same 
time, the scale of variation for Ax{y, z) along y off sample edges is macroscopic (see expressions (11))- Therefore, 
expansion (6) is also possible off the edges. Near the edges, however, Aa;(y, z) 7^ over entire thickness, and 
the scale of A{y,z) variation (see fig. 6b, d) in this region is A (in the y direction). Hence, expansion (6) is not 
valid in the limits y — b<y'<y + b near the edges (|y| —^ 1) of thick superconducting film. 
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3 The conditions for vortex entry in superconducting films 

Using expressions (8,11) it is possible to estimate the edge barrier suppression field s or the first- vortex entry 
field into superconducting films. 



3.1 Thin film 

The vortices may enter into a thin film in the Meissner state provided the condition |A(±1)| = Acrit is met; 
here Acrit ~ 'l'o/27r^ |l^ (<i>o is the quantum of a magnetic flux, is the coherence length). The resultant 
expression for Hg is 

with /3 being the same as in the expression (8). Dependence (12) in the limit W ^ Ae// and W ^ Ag// 
coincides to a factor of order one with the expression for the Meissner state breakdown field obtained in the 
limiting cases of wide and narrow thin films in ||l|. From (10,12) we can easily find the value of the magnetic 

field (or, rather, the thickness-averaged z-component) at a film edge H'f^^'^ , when vortices start penetrating in 

it. For example, for wide films ^ is 



To an accuracy of the factor of order one the above expression coincides with the field in the core of a Pearl 
vortex which is equal to $o/47r^Ae// [p^ . 



3.2 Thick film 

The main difference between thick and thin films is that the vector potential for the former largely depends 
on z. However, we should apparently anticipate that the conditions of vortex entry here will be qualitatively 
similar to those for thin films. Indeed, as shown in ||l^, |l^, after the vector potential has reached its critical 
value at a superconductor edge in the Meissner state (in the mixed state it is the gauge-invariant potential 
n = <I>oV(^/27r — A that should reach a critical value), the order parameter "if = |\E'|e"^ is strongly suppressed 
and vortex formation begins. The above papers dealt with bulk superconductors and thin-film samples, which, 
due to the symmetry of the problem or problem geometry, could be assumed homogeneous along the z-axis. 

It should be expected that in thick films the order parameter will be suppressed in the regions where the 
vector potential Ax{y, z) reaches its critical value. 

First, the condition |Aa;(?/, z)\ — Acrit will be satisfied at the side edges (±1, ±1) of a superconductor (as the 
magnetic field is increasing from zero). It means that the order parameter will be suppressed in these points 
first. With a further increase of the magnetic field this situation may develop by two scenarios: 

1. Suppression of the order parameter results in tilted vortices that start to form at the corners of a 
superconductor cross-section. When the vector potential at the equator reaches the critical value, two tilted 
vortices fragments (from top and from bottom) will join each other to form one rectilinear vortex. In the absence 
of pinning the latter is able to penetrate into the sample centre driven by the Lorentz force. Similar vortex 
entry scheme was discussed in 

2. In the course of further magnetic field increase, the order parameter becomes suppressed in the region of 
side edges. This causes areas with a suppressed order parameter to appear near the side edges, which would 
allow to regard a film cross-section as a rectangular with rounded-off edges. It should be emphasized that 
the geometrical sizes of sample remain unchanged in this situation, only its physical properties vary in the 
regions near side edges. This scenario allows to explain the physical mechanism behind the formation of the 
"geometrical rounds-ofF' near the corners of a rectangular cross-section sample, which were considered in ||]. 
However, unlike |3[ our approach is based on the assumption that vortices will start entering deep inside a 
superconductor when the condition |yla;(j/, 2;)| = Acrit is fulfilled at the equator. By that moment the effective 
"round-off" radius will reach a value of order d/2. 

Which of these two scenarios is practically feasible can be found out only through experiment. Theoretically 
this question can be answered by numerical solution of a problem on a vortex entry in a 3D sample within the 
nonstationary theory of Ginzburg-Landau. 
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The feature common for the above two schemes is, actually, the assumption that vortices enter deep a thick 
film when the condition |Aj;(l,0)| = Acrit is met. This allows to estimate the field of vortex entry inside a 
sample. Using Eq.(f f a) (regardless of possible variation due to the penetration of tilted vortices or the existence 
of areas with a suppressed order parameter), we now derive the expression for field Hs, which is equivalent (12) 
(with (3 = 2Xef f /ttW). This similarity is due to the fact that the Ax{il,0) value is defined practically by the 
same expression for both thin and thick films, provided parameter X^/dW is the same. 

By analogy with thin films, one can find a value of the field at a film edge when the vortices start to penetrate 
deep into a sample. Yet, as opposed to thin films, field H^'^^'^ largely depends on the z coordinate in the sample 
region. Estimations of the equator field yields 



which is practically equivalent to the H'^f^'^ (see fig. 7). Note also that H^'^{Hs) to the factor of order unit is 
equal to the thermodynamical field He, or the surface barrier suppression field for bulk superconductors in the 
absence of edge defects. 



3.3 The influence of surface defects and anisotropy 

The resultant values for field Hg are characteristics of isotropic superconductors with ideal surface. As was 
established in |jl^, |l4j, surface defects can considerably decrease the value of Hg. For example, in Q for the 
case K 3> 1 (k is the Ginzburg-Landau parameter) maximum suppression of the entry field g = Hg/ H^n (where 
Hf,n is the field of vortices entry in a superconductor with surface defects) was estimated as g ^ k/tt. Thus, 
for K, = 100 we will have g ~ 5.6. 

A strong influence on the value of Hg may be produced also by anisotropy or, rather, layered structure of 
superconductors. If the layers are not Josephson-coupled (or are weakly coupled), a superconductor should be 
regarded as a stack of superconducting layers. This geometry can be simulated, if we multiply the integrand 
in equation (3) (and, hence, in integral (6)) by the step periodic function z which is equal to one in the 
superconducting layer and is zero in the interlayer space. 

Let a period in such a structure be much smaller than A (which is practically always fulfilled for HTSC), 
a layer thickness be I and an interlayer separation be m. Then we can assume the distribution of A^iz) to 
be a smooth function z, similar to the dependence Ax{z) for a homogeneous superconductor. In this case, 
the integral in Eq.(3) for a layered superconductor will be [l + m)/l times smaller than that for an isotropic 
superconductor. In other words, we can replace this integral for a layered superconductor by the integral for 
an isotropic case by introducing an effective penetration depth A' = A-y/ {I + m)/l. In this way we immediately 
obtain the distributions of A{y), ^^^(y, 1), and also the values for Ax{y,z) at the side edges and the equator 
with allowance for the layered structure of superconductor. 

One particular effect of the anisotropy is that the value of ^2,(1, 0) will be \/ (I + m)/l times larger than that 
for an isotropic superconductor, all other conditions being equal. This, in its turn, will lead to a {I -\- m)/l 
times smaller field of vortex entry in a superconductor. For example, at I = 3A and m = 12A, typical for 
BiBaCuO, one finds ^J{l + m)/l k. 2.2. 

Thus, the two above-mentioned factors, i.e., surface defects and layered structure of superconductors may 
cause a considerable (10-fold and more) decrease of the vortex entry field in layered superconductors with surface 
defects, as compared to the vortex entry field in isotropic perfect-surface superconducting films. 

Another conclusion following from the fact of a layered structure in such systems is that the scale of a 
local magnetic field penetration, for example, for thick films, will be A' = (/ + m)/l, and for thin films the 
parameter Ae// has to be changed for A^jj = A'^/d. At the same time, the thickness-averaged current density, 
unlike the thickness-averaged vector potential, will practically remain unchanged across the entire film width, 
except for the regions lying close to sample edges. This can be accounted for by the fact that the expression 
for current density j{y) = ~cA{y) / AttX^ will include A only at a superconductor edge (see (8)). Likewise, all 
other quantities that obviously do not include A (for example, a degree of the magnetic field concentration at 
the equator of thick film) will remain the same. 
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4 The structure of a mixed state 



Let us now discuss the parameters of a mixed state arising in a superconducting film in fields larger than Hs- 
Here we neglect the presence of bulk pinning, which is justified for soft superconductors. This problem was 
studied earlier in ||, |[ ||, |l5| . In 0, ^ the authors considered a case of narrow thin films Wd/\^ ^ 1, 
1^ deals with wide thin films Wd/X^ ^ 1, and H, |^ is a study on thick films that formally obey the condition 
Wd/ » 1. We will analyse a general case to show that it embraces either of the above two situations. Besides, 
the resultant distribution of current density will be finite across the entire film width, as opposed to the results 
of the above cited works. 

Consider a superconducting film in a mixed state, placed in a perpendicular magnetic field. The current 
density and the vector potential in the London limit are related as j = —(A — $oV(/?/27r)c/47rA^, where cp is 
the order parameter phase. The Maxwell equation will have the form (1), in which A is now a 3D Laplacian 
operator. Using the Green function for A, we can write this equation in an integral form: 



A(r) = AO(r) + ^ / / / y^^dx' dy' dz' . (13) 

We now subtract function V(/3(r) from the left- and right-hand parts of Eq.(13) and take curl from these 
parts (using the property V x V</j(r) = 27r(5(r — r'), where r' — {x' , y') are the vortex coordinates). Next, we do 
the averaging over coordinates x, y on scales much larger than an intervortex distance. After these operations, 
the distribution of a current density becomes uniform along the x-axis and we can perform integration over x' 
in (13). Then we average the obtained equation over a film thickness and use the results of the integral (6) 
calculations. This will yields an equation for the sheet current density i{y) = f^^^^^j^iy, z)dz, 

+ - f ^dy' = -H^+niy)^o, (U) 
cW dy c J_^y- y' 

where n{y) is the density of vortices, the distance being measured in units of W/2. For the first time this 
equation was derived in Q| for thin films c? <C A. Just like Eq.(7), (14) is valid at |y| < 1 — \b\ for thick films, 
and across the width for thin films (excluding an extremely narrow near-edge region of size d ^ A). Besides, it 
should be expected by analogy with the Meissner state that Eq.(14) for thick films will also be valid directly at 
a sample edge. 

We analysed Eq.(14) numerically for different values of parameter W^/Ae//, using the condition that current 
density is zero in the region where vortices exist, and takes finite value in vortex-free regions |^, ^, |^. Besides, 
we set a boundary condition :/(±l) = ij,, on a current density (in increasing magnetic field), which allows for 
an edge barrier. The value of the current density js of order of the Ginzburg-Landau depairing current density 
for ideal-surface superconductors 12 1. In result, we have obtained the approximation expression for i{y) 



3 0<|y|<a, 
^^^y)^, , ^^-(" + 1) sign(y) a<|y|<l, (15) 



47r 



where 



2(2/2 - a^) 



-1, 



(l-« 



o-{Hoc) is the half- width of the vortex-filled region and parameter a is defined by expression (9) in which W 
has been replaced by W{1 — a). 

Expressions (15), being not derived, represents a rather useful interpolation for the distribution of the sheet 
current density for a film in a mixed state. We would like to emphasize again that in the thin-film case W/Xeff 
can be both smaller (narrow films) and larger (wide films) than unit, whereas for thick films this ratio is always 
much larger than 1. 
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Fig. 8 shows the dependence i{y) for different values of a magnetic field. The difference of approximation (15) 
from the numerical result does not exceed 4% in the vortex-free zone {a < y < 1). Note, that in the near-edge 
region and close to the boundary of the vortex-filled region deviation may come to about 9%. 

The dependence a(Hoo) (in increasing magnetic field) is to be found from the following expression: 

8 1 Xeff , 16 




For W/Xeff » 1 we have 
or 



while for W/\eff ^ 1 we have 

a{Hoa) = 1 — Hs/Hoa, 

for all values of Hryo ■ 

Using dependence (15), we can find distribution of the z-averaged magnetic field across a film width: 

HM^H^ + - t ^^dy'. (16) 

c J-iV- V 

Fig. 9 shows the dependences H z{y) for a film with W/\eff — 200 and a = 0.6 {Hoc — ^-iHs), obtained 
numerically and by means of expression (15,16). It is seen that these dependences practically coincide across 
the entire width of a film, except for the near-edge regions and the boundary of the vortex-filled zone. The 
difference in the magnetic field value aX y — 0.6 should be attributed to the inaccuracy of approximation (more 
precisely, its first derivative at the boundary of the vortex- filled area). Dependence Hz{y) shown in the insert 
to Fig. 9 was obtained theoretically in H, ^. It is seen that, as opposed to this analytical dependence, a 
non-zero magnetic field does exist outside the vortex-filled region also, and it is quite strong (> O.liJoo) for a 
film with the given parameters. Another distinguishing feature is the occurrence of a jump from zero to some 
finite value for the dependence n{y) = Hz{y)/^o at y ~ a. The reason of the vortex density discontinuity is 
explained by a non-zero magnetic field in the region (a, 1) and the condition of a magnetic field continuity at 

Using Eq.(15), we can estimate the dependence of /Hoc on the parameters of a film and an increasing 
external magnetic field (for thin films; see subsection 2.2) 




/3(1 - a) 
2^1 + ay2 



I -a' 



(17) 



It is nicely seen that in the limit a l(i?oo oo) the ratio ^ /Hoc tends to 1. 



Knowing the dependence i{y), we can find the magnetization curves of superconducting films for different 
values of W/Xeff. Fig. 10 illustrates the obtained results. One can see that with a increasing parameter W/Xeff 
the magnetization curves become similar to the dependence M(H) derived theoretically in ||^ for wide thin films. 
As W/Xeff decreases, the magnetization curves tend to the dependence which is valid for thin narrow films 
0,|l|l. Thus, even at W/Xgff = 1 the dependence M{H) for a narrow film and the M{H) obtained numerically 
by the use of Eq.(15) practically coincide. So, expression (15) allows to obtain magnetization curves for such 
superconductors at arbitrary ratio W/Xeff - Note that the magnetization curves in this case, i.e., at any value 
of parameter W/Xeff lie between two curves - 1 and 3 as shown in Fig. 10 (in dimensionless units). 
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5 Conclusion 



It is shown in the present paper that the Maxwell-London equation used so far only for thin films is also valid 
for samples of finite thickness. This equation is shown to define the distribution of a thickness-averaged vector 
potential and/or current density (in the mixed state case) across a sample width. For thin films the equation 
holds practically everywhere in a film, whereas in the thick film case its applicability is restricted only within a 
narrow bands near the edges — d/2 <\Y\ < W/2. 

An approximation expression is found, describing distribution of vector potential A (or current density j{y)) 
across the width of a film in the Meissner state, which applies to both thin and thick films. For thick films 
analytical expressions have been derived, defining the value of the vector potential (local current density) at the 
equator {Y = 1,Z = 0), side edges {Y = 1, Z = 1), and also on the top and bottom surfaces (Y, Z = ±d/2) of 
a sample. Besides, analytical approximation expressions have been found for the magnetic field at the equator 
and for the thickness-averaged edge magnetic field. 

The vector potential distribution data were used to evaluate the field of the the first vortex entry (barrier 
suupresion field) for superconductors of such geometry. It's described by an universal expression (12) valid for 
both thin and thick films. It is shown that besides surface defects the layered structure of superconductor may 
result in a significant (up to 100%) suppression of the vortex entry field. Thus, mutual influence both surface 
defects and layered structure may lead to suppression of Hs by factor ten (and even greater). 

The study of a mixed state yields an interpolation expression for distribution of a sheet current density i{y) 
in superconducting films without bulk pinning. This result allows to for the first time estimate the dependence 
on the magnetic field concentration 7 = H ^ ^ / H^o on the parameters of superconductor and external magnetic 
field ifoo- Besides, these data were used to calculate the magnetization curves for film superconductors at any 
values of parameter Wd/}?. 
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Figure captions 
fig.i 

Geometry of the problem. 
fig.2 

Distribution of averaged vector potential, for different ratios M^/Ag//: curves 1-5 for W^/Ae// — 1,5,10,50,200, 
respectively. Dotted lines are numerical results, solid lines are approximation (8). 
fig-3 

Dependence of the parameter a on W/Xeff. Circles are numerical results, solid curve 1 is the approximation 

(9). 

fig.4 _ 

Distribution of Hz inside the film for W/Xeff = 200. Dots are numerical result, solid line is expression 
obtained from approximation (8). Insert shows dependence H'I'^^'^ {W/ Xef j)'- circles are numerical result, dotted 
line is the expression (10), solid line is the fitting function 1 + 0.66^yW/Xeff■ 

fig.5 

Distribution of averaged vector potential at thick film {d — lOA) for various widths: W = 50A {1),W — lOOA 
(2), W = 200A (3). Dots are numerical result, solid lines are expression (8) with a — 0.25, (3 = 2}? jiidW . 
fig.6 

Contour lines of the intensity of magnetic field (a,c) and current density (b,d) of thick (W^ = lOOA, d ~ lOA) 
film in applied perpendicular magnetic field. The step for magnetic field is 0.41i7oo, for current density is 
0.1j(l,l). Maximum values of magnetic field (i? — 4.1iJoo) and current density [j — 1) are reached at the 
corners of the film. 

fig.7 _ 

Dependences H'^^^'^ (circles) and Hz{l, 0) (dots) on the parameter ^/wJd. Solid line 1 is the fitting function 
1/3 + 1.03y/W/d, dotted line 2 is the fitting function 1/3 + 0.92^/W/d. Insert shows the dependences of 

Hz (circles), iJ^(l, 0)(dots) and /^^(l, 1) (stars) are shown on the film thickness for W/d — 5. 
fig.8 

Distribution of the sheet current i{y) for film with parameter W/Xeff = 200 and for different values of a: 
0.0 (1), 0.4 (2), 0.8 (3). Dotted fines are numerical results, solid lines are expression (15). 
fig.9 

Distribution of the averaged z-component of magnetic field for W/Xeff = 200 and a — 0.6. Sohd fine is 
obtained with help of expression (15,16), dotted line is numerical result. Insert shows detailed distribution of 
the field; dashed line is the function H \faF~^]p / \p\~^]p from ^. 

fig.lO 

Magnetization curves of superconducting films for different ratio W/Xeff'. curve 1 for W/Xeff — oo, curve 
2 for W/Xeff = 200, curve 3 for W/Xeff — 1. Curve 3 practically coincides with the magnetization curve for 
narrow W ^ Xef / films ||T^ . 
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fig. 1 D.Yu.Vodolazov and I.L.Maksimov 




fig. 2 D.Yu.Vodolazov and I.L.Maksimov 
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fig. 3 D.Yu.Vodolazov and I.L.Maksimov 




fig. 4 D.Yu.Vodolazov and I.L.Maksimov 




fig. 6 D.Yu.Vodolazov and I.L.Maksimov 




fig. 7 D.Yu.Vodolazov and I.L.Maksimov 




fig. 8 D.Yu.Vodolazov and I.L.Maksimov 




fig. 9 D.Yu.Vodolazov and I.L.Maksimov 




fig. 10 D.Yu.Vodolazov and I.L.Maksimov 



